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ABSTRACT 

We describe a new mm-wave molecular-line mapping survey of the southern Galactic Plane and its 
first data releases. The Three-mm Ultimate Mopra Milky Way Survey (ThrUMMS) maps a 60° x 2° 
sector of our Galaxy’s fourth quadrant, using a combination of fast mapping techniques with the Mopra 
radio telescope, simultaneously in the J=1^0 lines of ^^CO, ^^GO, C^^O, and CN near 112 GHz at 
^arcminute and ^0.3kms“^ resolution, with ^2 K channel”^ sensitivity for and ^1 K channel”^ 
for the other transitions. The calibrated data cubes from these observations are made available to the 
community after processing through our pipeline. Here, we describe the motivation for ThrUMMS, 
the development of new observing techniques for Mopra, and how these techniques were optimised to 
the objectives of the survey. We showcase some sample data products and describe the first science 
results on CO-isotopologue line ratios. These vary dramatically across the Galactic Plane, indicating 
a very wide range of optical depth and excitation conditions, from warm and translucent to cold and 
opaque. The population of cold clouds in particular have optical depths for easily exceeding 

100. We derive a new, nonlinear conversion law from integrated intensity to column density, 

which suggests that the molecular mass traced by CO in the Galactic disk may have been substantially 
underestimated. This further suggests that some global relationships in disk galaxies, such as star 
formation laws, may need to be recalibrated. The large ThrUMMS team is proceeding with several 
other science investigations. 

Subject headings: astrochemistry — galaxies: the Milky Way — ISM: kinematics and dynamics — 
ISM: molecules — radio lines: ISM — stars: formation 


1. INTRODUCTION 

The Milky Way is a large, massive, and luminous 
disk galaxy, and since we are immersed in it, it can 
provide detailed and localised keys to understanding 
larger-scale processes that operate in other galaxies 
throughout the observable universe. For example, our 
Galaxy uniquely enables parsec-scale opportunities to 
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cesses which are typically not resolved to better than 
^100 pc scales in even nearby galaxies. 
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Therefore, high-resolution wide-held surveys of the 
Milky Way have an important impact on several of these 
topics. For example, in studies of star formation, there 
is a critical threshold for physical resolution of about a 
parsec, since that is the relevant scale for the formation 
of individual clusters from dense molecular clumps, and 
thes e are likely the basi c unit of star formation activ¬ 
ity (Lada & Lada 2003) throughout the Galactic disk. 
For a typical distance of 3kpc, the clump/cluster scale 
is only reached with angular resolutions of an arcminute 
or better. 

At the largest scales, rece nt wide-held Galytic P lane 
surveys such as GLIMPSE (i Benjamin et al . ||2003 |) and 
its successors, ATLAS GAL (Schuller et al.JzOOOl), and 


Hi-GAL (Molinari et al. 
understanding of the IVIil 


20I0|) are revolutionising our 
iy Way, by enabling the global 


study of the cold molecular structures of the interstellar 
medium (ISM) and the star-forming activity within these 
structures, throughout the disk. Most, if not all, of these 
surveys have been made feasible by the development of 
wide-format cameras or multi-pixel bolometers, enabling 
the coverage of hundreds of square degrees of sky at an¬ 
gular resolutions of only a few arcseconds. While sen¬ 
sitively revealing the two-dimensional morphology and 
broadband spectral energy distributions (SEDs) of many 
objects, many of these surveys are in the continuum only; 
that is, they cannot provide kinematic nor line-of-sight 
distance information. 

Comprehensive spectral-line surveys with this level of 
wide-field coverage are rarer and much more limited: the 
main wide-field surveys of note are in HI (e.g., the SGPS 
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Figure 1. Columbia-CfA survey Wco image of the Milky Way’s 4th Quadrant between longitudes 359° and 299° ( jPame et al.| 

|200l|. The black-to-red colour scale is over 0-200Kkms“^. The red box shows the original planned coverage of ThrUMMS. As of 2014 
Dec, mapping of this region is ~80% complete, yielding ~7x higher resolution maps than the CfA survey (~50x higher areal resolution), 
and simultaneously in the 4 molecular lines ^^CO, ^^CO, C^^O, and CN. 


in the south; 
the Columbia-CfA survey: 


McClure-Griffiths et aL||2QQ5 ), ^^CO (e.g., 


Dame et al. 


2001 ) or i^CO 


(the GRS in the north; | Jackson et al.||2(J06|), and mor e 
recently H 2 (UWISH 2 in the south; Ihroorich et al.| 2011 |) . 
Fewer still have adequate spatial resolution across these 
wide scales. So while detailed surveys of portions of the 
Galactic ISM exist, each of these necessarily has its set 
of limitations. 

Thus, UWISH2 lacks a practical velocity resolution for 
comparison to the radio-frequency surveys, and focuses 
on only one spectral line, limiting its power as a physical 
diagnostic. The SGPS is wide-held, but the velocity and 
angular resolutions are too modest (lkms“^, 2 ') for de¬ 
tailed work involving correlation and comparison of the 
numerous and signihcantly higher resolution near- to far- 
IR datasets mentioned above. Further, the HI mapped 
by such surveys comprises a different phase of the ISM 
than that directly involved in star formation. The ven¬ 
erable Golumbia-GfA survey (Fig.[^ is at a resolution of 
only 8 ', or ^7pc for a typical “nearby” (i.e., distance 
< 3kpc) molecular cloud, and is only in one species. 
The Nanten survey is only slightly higher resolu¬ 

tion (4'), s till well above the cffi mp resolution threshold. 


I 2 OO 6 D is c 


The GRS (Jackson et al 
tion ( 0 ( 8 ), but only cover's 
and again i s in just one species (^^CO). GHaMP (Barnes 
et al. [ 2011 ) covers 120 deg^ of the fourth quadrant and 


at a reasonable resolu- 
of the hrst qu adrant. 


is multi-species and reasonable resolution (Offi), but the 
helds mapped are of moderate width (^ 1 ° or less) and 
focused on the denser molecular clumps. There are also 
a few other ^degree-wide molecular cloud maps in the 
fourth quadrant recently made at Mopra, but their areal 
coverage is substantially less than GHaMP’s. 

Thus, by various metrics, we have no complete, wide- 
held, multi-species inventory of the gaseous content of 
the ISM at an angular resolution approaching that (20" 
or better) of the continuum surveys. Fut ure large-scale 
Galac tic ISM surveys, such as GASKAP (|Dickey et al.| 
2Q13[ ) , will soon yield HI maps with resolutions of 20 " and 
O.lkrns”^ and high sensitivity, providing a signihcant 
improvement over the existing SGPS and its sister HI 
surveys. But to explore the molecular content of the 
Galaxy, we must still contend with some fundamental 
observational limitations. 

ThrUMMS’ genesis occurred in 2009 at a meeting 
of the Mopra community (the Millimetre Astronomy 
Legacy Team, or MALT), where a need for three large 
legacy projects was identihed, each covering the southern 
Galaxy in a different frequency range and having related. 


but somewhat different, science goals (see below). The 
110 GHz project (subsequently dubbed ThrUMMS by its 
participants, and also called MALT- 110 ) was designed to 
address the wide-held issues described above. 

In this paper we introduce the Three-mm Ultimate 
Mopra Milky Way Survey (ThrUMMS), and provide the 
background information on ThrUMMS for the commu¬ 
nity to make best use of ThrUMMS’ wide-held, molecular 
transition data products. In the next section we describe 
the survey design and the unique operational aspects of 
the data collection and reduction. In ^we describe the 
current Data Releases and the data products they con¬ 
tain. il presents our hrst science results for the global 
dataset, ^concludes with a brief summary. 


2 . SURVEY DESIGN 

2 . 1 . The Need for Faster Mapping Teehniques 

Even today, most radio telescopes have single-pixel de¬ 
tectors in their focal plane. Eor single-dish antennas op¬ 
erating over the ^ 1-100 GHz range, angular resolutions 
are typically limited to about an arcminute. Without 
heterodyne focal-plane detector arrays (which are still 
rare and have most recently tended to focus on narrow- 
held astrochemical studies) to speed up the mapping, 
even this low resolution takes an extraordinary time to 
completely map large areas of sky at reasonable sensi¬ 
tivity. The conventional wisdom for some time has been 
that the dedicated, small-telescope, GfA/Nanten model 
was the only practical way to produce useful wide-held 
molecular maps. 

Motivated by this conundrum, we devised two new 
techni^es specihcally optimised for the 22 m diameter 
Mopr an antenna, to dramatically improve the telescope’s 
speed for wide-held mapping. We also implemented a 
third technique (pioneered in the GfA survey) that im¬ 
proves the image hdelity for all mapping projects. A 
detailed exposition of th ese techniques has been given by 
Barnes & Dame (2009), but for completeness we high- 
light the essential features here (summarised in Table [h. 


2.1.1. Fast Mapping with Mopra 

The hrst of these techniques, “East Mapping” (EM), 
takes advantage of the fact that, especially for the ^^GO 
line but also for ^^GO, the S/N performance enabled 


® Operation of the Mopra radio telescope is made possible by 
funding from the National Astronomical Observatory of Japan, the 
University of New South Wales, the University of Adelaide, and the 
Commonwealth of Australia through CSIRO. 
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Galactic Longitude 


Figure 2. RMS noise maps for a combination of BSM and AM 
(above, from Dame et al., in prep .) and normal OT F mapping 
with a fixed (right, from|Barnes et al.||2011|), each over 

a large number of OTF fields (26 & 16 resp.j. i ne ettective beam 
sizes for the different mapping techniques are shown in the lower- 
left and lower-right corners, resp. Both maps are shown on the 
same brightness scale to make the comparison more straightfor¬ 
ward. For the AM/BSM map, there is only one modal peak in the 
rms at 0.65 K, due to the sky-active setting of ujsf. For the normal 
OTF map, the modal rms is 0.75 K, but there is a much larger 
variation in rms than in the AM/BSM map due to the rms’ strong 
elevation dependence. There are also blended secondary peaks in 
the histogram (shown in Fig.l^ at 0.85 K due to poor weather for 
one OTF field, and at 0.65 Iv due to excess coverage for another 
field. 


by a 22 m diameter dish is very high, and for survey 
mapping of such bright lines, sensitivity requirements 
can be significantly relaxed compared to, for example, 
those needed for the common “dense gas tracers” such 
as HCN or CS. For the latter, at Mopra the sampling 
interval of 2.0 s coupled with Nyquist-sampled on-the-fly 
(OTF) mapping produces maps with typical rms noise 
of 0.3Kch“^, which is well-suited to mapping of such 
lines with characteristic brightness ^1-3 K. However, in 
the radiometer equation, the rms noise arms for 

integration time t. So for example, a factor of 3 worse 
sensitivity, i.e., a ^ lKch“^, for a very bright line (5- 
20 K or more) dramatically reduces the integration time 
necessary to achieve a still excellent S/N, by (in this ex¬ 
ample) a factor of 3^ ^ 10. In other words, with suitable 
modifications to the observing hardware and software, we 
can significantly reduce the sampling interval in a given 
observing mode, such as OTF mapping, and allow maps 
of bright lines to be completed perhaps an order of mag¬ 
nitude faster than with regular “slow” mapping. In prac¬ 
tice, the FM sampling interval was fixed to be about 8x 
faster than the regular sampling interval, namely 0.256 s 
instead of 2.000 s. 


2.1.2. Beam Sampled Mapping with Mopra 

The second technique also takes advantage of a par¬ 
ticular feature of Galactic emission: it is spatially 

very widespread. Again, for a survey mapping project, 
Nyquist sampling of the emission (especially at 

the more modest sensitivity of FM) can be considered 
slightly superfluous. Conveniently, the OTF mapping 
software at Mopra has an easily-set parameter which con¬ 
trols the degree to which the spatial sampling matches 
(or doesn’t match) the Nyquist-sampling criterion, called 
the “Nyquist frequency,” un- Thus, when (the 
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observing frequency), sampling during regular OTF 
mapping is set by default to be ^10% finer than the for¬ 
mal Nyquist spacing at Uq. For ThrUMMS, however, we 
set v^—Voj2^ so that OTF mapping in this case is per¬ 
formed with the telescope scanning the sky twice as fast 
as the default. This effectively means that ThrUMMS 
maps are sampled at ^half the Nyquist rate (i.e., in a 
spatially coarser way), in both the scanning direction and 
between scan rows. With this change, OTF mapping ob¬ 
tains effectively beam-spaced (33" at 115 GHz) samples, 
a method we have dubbed “Beam Sampled Mapping” 
(BSM). This follows a long tradition of “complete” but 
beam-spaced mapping from the early days of ra dio as 


tronomy, including parts of the GfA survey itself (Dame 


Because this is a rather novel mapping approach, we 
describe here some aspects to BSM that, while not mean¬ 
ingfully affecting the data quality, are worth bearing in 
mind for some applications, e.g., those that depend on 
accurate structural information at the smallest effective 
resolution of the data. 

During normal OTF mapping, the antenna is contin¬ 
uously slewed across a raster line in the eventual OTF 
map. During the slew, data are accumulated almost con- 
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Figure 3. RMS noise histograms (log scale for y, same linear brightness scale for x) for the AM/BSM (left) and normal OTF (right) maps 
shown in Fig.|^ The narrower spread of system noise in the left panel is due to the noise-homogenisation effects of the Active Mapping 
technique. The slight difference in normalisation of these plots is explained in the caption to Fig.m 


Table 1 

Mopra mapping techniques 


Mapping 

Mode 

Sampling 

Interval 

Finest 

Resolution 

Atmospheric Stability 
Accommodated 

Relative 

Mapping Speed^ 

Relative 

RMS Noise^ 

regular 

2.000 s 

33" 

cannot adjust 

1 

1 (var.) 

FM 

0.256 s 

33" 

cannot adjust 

8 

\/8 (var.) 

BSM 

any 

66" 

cannot adjust 

4 

1 (var.) 

AM 

any 

any 

>lhr 

<1 

1 

ThrUMMS 

0.256 s 

66" 

>lhr 

<32 

U8 


^ These factors are multiplicative: for ThrUMMS, all of FM+BSM+AM are used. 


tinuously in the MOPS digital filterbank, and interrupted 
very briefly each sampling interval (normally 2.0 s) to 
dump out the contents of the spectrometer to mass stor¬ 
age. The speed of the slew is tuned to the sampling inter¬ 
val so that, for normal OTF mapping, Nyquist sampling 
on the sky (at the given observing frequency) is achieved 
by default (i.e., z^Ar=^o)- 

During BSM+FM, even though the sampling interval 
is now about 8 times faster than for regular mapping, 
the telescope is typically slewing across the map 16 times 
faster, because z^Ar=|z^o- Thus, the data dumps from the 
spectrometer occur about once per beam in the raster 
direction. Because the integration is continuous over 
this time, we are effectively forming a properly-smoothed 
average spectrum along the raster direction at ^beam- 
spaced positions, and apart from this spatial convolution, 
there is no loss of information (such as zero-spacing data) 
on source structure on these scales. 

Across the raster lines, however, the situation is differ¬ 
ent. Here the selection of means that the raster 

lines are separated by one intrinsic Mopra beamsize, and 
we are therefore undersampling the sky emission in this 
direction. In other words, we could be missing some zero¬ 
spacing information in this dimension. However, Mo¬ 
pra’s intrinsic angular resolution of 0(6 at 3 mm is well 
below the resolution threshold for cluster-formation (see 


P , so even with beam-spaced sampling, this issue is not 
particular concern for ThrUMMS science. For exam¬ 
ple, in some planned science investigations, ThrUMMS 
maps are intended to be starting points for more de¬ 
tailed follow-up studies. More importantly, we allow for 
this feature of the mapping in our data pipeline, produc¬ 
ing reliable ma ps of the emission at the quoted effective 
beam size (see ^2.5). 

With BSM, since we do not change any of the parame¬ 
ters in the radiometer equation, we pay no noise penalty 
per se for such maps. Obviously, however, the maps so 
generated are at twice Mopra’s inherent angular resolu¬ 
tion. But mapping a given area now proceeds four times 
faster than before, which for ThrUMMS we consider a 
very valuable trade-off. In summary, the net effect of 
combining the FM and BSM techniques results in maps 
of sensitivity ^I K ch“^ at an angular resolution of ^1(2, 
but made at an areal mapping speed roughly 32 x faster 
than Nyquist-sampled “slow” OTF mapping. This is a 
dramatic improvement over (e.g.) the CfA or Nanten 
maps, and ideal for our purposes. We call this combina¬ 
tion of techniques “Very Fast Mapping,” or VFM. 


2.1.3. Active Mapping with Mopra 

The third technique, new to Mopra but used at some 
other facilities, is a sky-related modification of BSM, 
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Table 2 

ThrUMMS 3mm band spectral lines 


Species 

Transition 

Frequency 

(GHz) 

Efficiencies^ 

m rjc 

Resolution® 

(kms-i) 

RMS'^ 

(K) 

Science Topics 

12CO 

J=1^0 

115.271 

0.42 

0.55 

0.33 

1.3 

largest scales, GMC envelopes, temperatures, outflows 

CN 

J=1^0 

113.5 (5 hf)d 

0.42 

0.56 

0.33 

0.9 

dense clumps, Zeeman candidates 


J=1^0 

110.201 

0.43 

0.57 

0.34 

0.7 

PDRs, fractionation, optical depth, velocity dispersion 

Ci®0 

J=1^0 

109.782 

0.43 

0.57 

0.34 

0.7 

column density, depletion, dense gas finder 


^ Following |Kutner & Ulich| (|1981|), the beam efficiency 77 ^ is appropriate for obtaining T^b of compact (~beam-sized) sources, while 
the coupling etticiency 77 c is appropriate for computing T^b of more extended (several arcmin) sources. 

^ This is for a 4-channel binning: the intrinsic velocity resolution is 0.085kms“^. 

This is for a 4-channel binning: the intrinsic rms noise per full-resolution channel is twice as large. 

^ This line has 9 hyperfine components, of which 5 are recorded in the IF zoom band. 


but which cau also be applied to regular “slow” OTF 
uiappiug. For uiost uioderu uiui-wave radio telescope 
systeuis, the systeui teuiperature Tgys is douiiuated by 
the coutributiou froui the atuiosphere, Tatm- Apart froui 
weather-related factors such as huuiidity, this of course 
depeuds ou the airuiass at the elevatiou e of the ob- 
servatiou above the horizou, Tatm cx: cosec(e). It turus 
out that we cau easily (aud exactly) accouiuiodate the 
predictable, elevatiou-depeudeut couipoueut of Tgyg by 
scaliug un to Tgyg(e)/Tgyg(90°), siuce the slightly slower 
uiappiug at a larger will iucrease the iutegratiou tiuie 
iu exact proportiou to that required to adjust for the 
higher Tgyg(e) couipared to Tgyg(90°). Thus, 


z/iv(e) 


*^iv(90°) 


/ Tsysjs) \ 
\T,ys{90°)J ’ 


( 1 ) 


with the uiiuiuiuui z^Ar(90°) choseu to be 55 GHz. We call 
this techuique “sky-active uiappiug”, or AM. 

Iu practice, the uiiuiuiuui Tgyg(90°) for the scaliug of 
un was choseu to be a couservative approxiiuate global 
uiiuiuiuui iu the best couditious (400 K at 115 GHz iu 
our case), rather thau estiiuated separately for each day 
at e=90°. Iu that case, slowly-varyiug (^hourloug) aud 
day-to-day weather chauges could also be accouiuiodated 
by this scheuie just as easily. Ouly iu the uiost vari¬ 
able weather is the AM luethod susceptible to bad spec¬ 
tral baseliues, produciug “weather stripes” iu the luaps, 
but this is true for auy OTF uiappiug project, aud the 
weather iu such cases is, ou average, uot good euough 
to observe at all. This sky-active techuique at least 
eveus out luuch of the weather variatiou, aud usually 
the scheuie works quite well: au illustrative exauiple of 
the couibiuatiou of two of these techuiques (AM+BSW 
usiug Mopra’s OTF systeui is showu iu Figures 2 aud[^ 
aud is couipared with a siiuilar luosaic luade with regular 
OTF uiappiug. 


2.2. Spectroscopy 

Mopra gives us auother iuiportaut advautage over 
uiappiug surveys luade elsewhere: the MOPS spectroui- 
eter. This further couibiuatiou luakes ThrUMMS rather 
uuique couipared to previous luolecular surveys of the 
Milky Way. 

For “slow” OTF uiappiug, MOPS siiuultaueously 
records up to 16 high spectral resolutiou (^0.1kuis“^) 
IF “zooiu luodes” of width 137 MHz each (^400kuis“^ 
at 3uiui), which together cau seiui-arbitrarily cover au 


8 GHz wide portiou of the 3 luiu baud. Especially iu the 
liue-rich 90 GHz rauge, Mopra cau therefore luap 16 or 
uiore spectral hues extreuiely efficieutly, as iu CHaMP 
( Barues et al.|2Qll ) aud other surveys of deuse luolecular 
cores aud clumps. Duriug FM, however, the data traus- 
fer rate froui the spectrouieter hardware to luass storage 
is luaxiiuised, aud iuiposes a liiuit of 4 IF zoouis ou the 
data-takiug. This liiuitatiou, however, does uot affect 
the scieuce goals of ThrUMMS, siuce there are ouly 4 
trausitious bright euough (at the seusitivity of FM) with 
a practical astrophysical applicatiou iu the 107-115 GHz 
rauge. These trausitious are listed iu Table We show 
iu Figure [^souie sauiple spectra to illustrate the spectral 
data quality duriug VFM+AM. 


Glon: 333.3104° 
Glat: -0,3591322° 



Velocity (km/s) 

Figure 4. Sample spectra from ThrUMMS data cubes. After 
a 4-channel binning (changing the intrinsic 0.085kms“^ channel 
width to 0.34kms“^), the per channel rms noise in each spectrum 
is ~1.3K (^^CO, top) and ~0.7K (^^CO, bottom). Note also that 
we typically have good baselines by choosing OFF positions at 
similar elevations to the observed maps. 
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2011 2012 2013 2014 


Figure 5. Calibration data from all f ive ThrUMMS obs erving 
“seasons” to date, on the source BYF99a (|Barnes et al.|2011|). Each 
point shows a fit to the emission in this source tor tJie calibration 
scan at the beginning of each day’s observing, on an ordinate scale 
normalised to the mean of all measurements. Shown in grey is 
the mean zb standard deviation of all these fits; the normalised 
standard deviation is 11%. 


2.3. Coverage of the Galaetie Plane 

In 2012 we completed (except as noted below) all map¬ 
ping for the equatorial degree-wide strip (i.e., \b\ < 0°.5) 
for two-thirds (360° > I > 300°) of the Fourth Quadrant 
(hereafter 4Q). In 2013 we began expanding the coverage 
to | 6 | < 1 ° over the same longitude range, and anticipate 
that this will be completed in 2015. While optimising 
the data-taking for efficiency and low noise, the filling of 
our 60° X 2° target area during the observing is somewhat 
random, but large contiguous portions of the 4Q are now 
available. A small number of fields are still affected by 
poor weather or equipment problems, such that they will 
require re-observation to maintain uniform sensitivity; 
other fields were not finished in one session. Despite this 
patchwork progress, there are no calibration “edges” be¬ 
tween fields observed at different times or seasons, since 


the calibration is robust (see ^2.4| and ^2.5). Moreover, 
because of our variable-approach to noise manage¬ 
ment, the rms varies only modestly between fields, with 
very little (if any) effect of such blended data being ap¬ 
parent in the merged data cubes. 


2.4. Observing Proeedures 

All ThrUMMS maps are made with the VFM OTF 
raster direction in longitude (i.e., raster lines are at con¬ 
stant latitude). The absolute pointing accuracy is mon¬ 
itored and corrected to within 6 " using various 86 GHz 
SiO maser sources (R Cy , IRSV1540, AH Sco, VX Sgr; 


Indermuehle et al. 2013) between each map. Pointing 
variations rarely exceed 20" at the end of each OTF 
map. (With the 24" pixel size in each data cube, such 
pointing corrections, although slightly larger than those 
for smaller OTF maps at Mopra, are inconsequential 
here.) At the beginning of each daily observing session, 
a calibration scan is taken of the bright molecular clump 
BYF99 in the 77 Carinae Nebula complex (Fig.l^, and 
the Tgys measurements made to start the sc^ng for 
subsequent VFM. OF F positions fo r each VFM field are 
taken from the list of Dame et al. (2001). Under good 
conditions each 0°.5 field takes about 2" of clock time. 


or somewhat longer as the Tgys rises. Fields are priori¬ 
tised for mapping so that, whenever possible, they start 
at an hour angle (HA) near - 1 ^, and finish near HA = 
+ 1 ^, minimising the air mass and maximising the map¬ 
ping efficiency. This procedure is controlled by a robotic 
telescope observing script, implemented in the control 
software since 2010 , and requiring only occasional over¬ 
sight functions from our observing team. 


2.5. Data Reduetion and Proeessing 

We use an au gmented version of t he Livedata/Grid- 
ZILLA package (|Barnes et al.||2001|), the usual standard 
for ATNF single-dish data reduction. This package pro¬ 
vides a convenient environment for basic spectral calibra¬ 
tion, baselining, separation of different IFs and spectral 
line data, gridding and convolution, and the formation 
of standard FITS (/, 6, V) cubes. We describe our proce¬ 
dures next, including custom improvements to the pack¬ 
age. 

In Livedata, raw data for each observed field are 
bandpass-corrected and scaled by the measured Tgys to a 
relative antenna temperature scale T^, using the “ refer - 
ence” method (with OFF positions as described in ^2.4). 
A linear baseline from a robust fit to signal-free channels 
(not including 100 channels at each end of each 4096- 
channel spectrum) is also subtracted, the frequency axis 


G328-340 



Figure 6. Integrated intensities from DR3 maps cover¬ 

ing 340°>/>328°and integrated over -55 <W SR <“40 kms“^, com¬ 
pared with t he same measurem ents from an equivalent area in the 
CfA survey (|Dame et al.|2001|). The DR3 data were convolved to 
the CfA beam ot 8T to torm tnis pixel-by-pixel comparison. (This 
is the same area and velocity range as displayed in the 2nd and 3rd 
panels of Fig.[^) Points from pixels at the edge of the ThrUMMS 
comparison area, where the convolution to the larger CfA beam 
becomes inaccurate, are excluded from this plot. The dashed line 
indicates a 1:1 equivalence; the average result over this and sev¬ 
eral other areas and velocity ranges is given by the red line at 
90zb5%. This indicates that the ThrUMMS data are on a bright¬ 
ness scale slightly below that of the CfA data, but that the relative 
consistency between the two surveys is excellent, vindicating our 
quality-control procedures. 
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Figure 7. Comparison images of a 12° section of the Milky Way’s 4th quadrant (4Q), as labelled. At 1^2, ThrUMMS’ high areal resolution 
allows us to explore the detailed physics of individual clouds and star-forming regions in much greater detail than before. Simultaneously, 
we retain the ability to address global Galactic structure topics with a wide-held, uniform, sensitive, multi-transition survey. In these 
pipeline-reduced ThrUMMS Data Release 3 (DR3) images, the rms noise is ~2-3Kkms“^, and we have lar gely removed the mapping 
artifacts present in DR2 and earlier data. Here we compare DR3 maps with equivalent Columbia-CfA survey ( [Dame et al. |2001| l data for 
two different velocity integrations across the 12° area. The upper panels (velocity range -55 to -40 kms“^) are compared numerically in 
Fig.|^ The higher resolution of ThrUMMS maps compared to the CfA maps allows easy separation of individual continuum features that 
reside at different velocities (and hence at different distances). 


is rescaled for all topocentric to Vi^sR corrections, and 
the result for each MOPS IF-zoom spectrum is written 
to a separate SDFITS file. 

Most of our custom augmentations to this pipeline are 
effected between the LivedATA and Gridzilla stages. 
First, absolute flux calibration is performed on the SD¬ 
FITS files produced by Lived ATA, as illustrated in Fig¬ 
ure Next, automated removal of known observing or 
processing artifacts is performed. For example, in DR2 
and earlier data releases, some artifacts (failed bandpass 
corrections for single 0.256-sec samples, bad calibrations 
due to observing malfunctions, or bad-weather “stripes” 
where baseline subtraction becomes less reliable) were 
not corrected, and were manifest in the data products. 


While the fraction of data so affected is small, the ef¬ 
fects on the cubes and maps can be quite large, since the 
individual errors can be of a large magnitude. Left uncor¬ 
rected, such data produce noticeable deleterious effects 
in the maps, such as “bright” or “dark” spots, strong 
“edges,” warped baselines, and so on. However, these er¬ 
rors were easy to identify and then either flag or correct. 
Thus in DR3, we have completely removed the single¬ 
pixel artifacts, systematic calibration errors, and some 
of the weather stripes, by new, custom, semi-automated 
routines written as a combination of IDL and Unix c-shell 
scripts. These additional calibration and editing steps 
improved the overall data quality to the point where the 
resulting noise levels were always consistent with theo- 
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Figure 8. Comparison of simultaneously-mapped ThrUMMS species on the scale, as indicated, integrated over -100 km s ^ < Vlsr < 
-h50kms“^. The contours in each image are the same, of at levels of 1, 2, ..., 7, and 8 times 13.6Kkms“^ 20cr. These DR3 maps 

cover half the area shown in FigureNote the different line ratios of various peaks, and that both and CN are easily detectable in 

several places. The and emission is quite filamentary, similar to Hi-GAL images. 


retical expectations, given the observing conditions. 

From this point the data processing resumes with the 
normal operation of Gridzilla, but with mapping pa¬ 
rameters that are guided by the scientific motivation of 
defining the Galactic CO and CN emission to arcminute 
and sub-kms“^ resolutions, given the unique features of 
the BSM technique. During the observations, one inte¬ 
gration cycle occurs approximately every half-arcminute 
interval, in both the scanning and scan-step direction 
(i.e., I then b). This is comparable to the native resolu¬ 
tion of the Mopra tele scope at 115 GHz (approximately 
33"; |Ladd et al.|p0^, SO our mapping strategy under- 
samples that scale. To form images with an effective 
resolution that is well-sampled, we define gridding pa¬ 
rameters in Gridzilla that are ^twice as large as those 
for normal OTF maps. Thus, the map pixels are of 
size 24" and the gaussian convolving beam is 68 ", trun¬ 
cated within a 36" gridding kernel, giving an effective 
resolution of 72" in the maps, with a typical contribu¬ 
tion of 40X0.256-sec samples per pixel to the final spec¬ 
tra. Therefore, the effective integration time is about 
lOspixeD^, and the output from Gridzilla preserves 
the intrinsic velocity resolution of MOPS, as listed in the 
footnote to Tablewith the typical noise level also given 
there. The final spatial sampling compares well with the 
typical spacing of ^20"-25" 1 pixel) between raster 

lines in our VFM OTF maps. Our mapping and data 


reduction procedures therefore return a product that is 
better than Nyquist-sampled at 72" resolution, and cover 
an area of the Galaxy not feasible otherwise. 

Table also summarises some parameters of the data 
cubes for the 4 spectral lines we are mapping. The per 
channel rms noise values drop with frequency since these 
lines are at the edge of the atmospheric 3mm band; thus, 
the Tgys for ^^CO and C^^O are almost a factor of 2 less 
than for the ^^CO data. Each species’ cube is produced 
on the scale, and must be converted to T*=Tmb for 
physical analyses. Table 1^ gives both the main beam ( 775 ) 
and inner error beam ( 77 ^ efficiencies we recom mend for 
eac h spec tral line, based on the measurements of |Ladd ^ 
ah] (| 2QQ5|). These efficienc ies are respectively appropri¬ 
ate ( [Kutner fc Ulich|[l981 ) for compact sources of ^one 
intrinsic Mopra beam or less, and for slightly more ex¬ 
tended (a ^few arcmin) sources. For most emission seen 
in the ThrUMMS maps, rjc will be more often used, and 
the DR cubes are all scaled to T* by r]c (although see the 
discussion in ^ 2.6 about refinements to this scaling). 


The next step in our augmented processing pipeline 
after the Gridzilla stage is to calculate moment maps 
for the line emission. We describe these methods in 83.31 


2.6. Verification 

The final step in the data reduction pipeline is to check 
that we have constructed maps that are at least consis- 
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Figure 9. Sample downstream products from the 
ThrUMMS DR3 cubes, where a smooth-and-mask 
(SAM) algorithm has been applied using a custom 
IDL-Miriad shellscript to calculate moment maps, 
(top-left) intensity, integrated over all un¬ 

masked voxels from -150 to +50 kms“^; contours 
are at 10, 30, and 50 K kms“^; (bottom-left) 
velocity field, with contours as in the intensity map 
and evaluated over the same voxels and velocity 
range; and (right) longitude-velocity diagram, 

integrated over all unmasked voxels in \b\ < 1°. Al¬ 
though artifacts due to poor weather are still visible 
in these data, these too will be removed in future 
DRs. Note in particular that beam-sized noise mod¬ 
ulation in the velocity field maps, common in such 
higher-moment calculations, are virtually eliminated 
in these DR3 products. This is due to the superior 
filtering properties of the SAM method. The effec¬ 
tive Mopra beam is shown in the bottom-right corner 
of each of the intensity and velocity maps. 


Galactic Longitude 


tent with p revious surveys, e.g., the CfA data of |Dauie 
et al. (2001). As can be seen in Figure the ThrUMMS 


data do indeed compare well: averaging over a variety of 
regions, we derive an overall scale factor between the two 
surveys of T; (ThrUMMS) = (0.90U0.05) x T;(CfA), 
which is within the overall calibration uncertainties given 
in FigureWe consider this result (close to unity, with a 
small variance) a strong verification of the quality-control 
procedures described above. 

Even so, the conversion factor is slightly but system¬ 
atically below 1 . This may be due to two effects. One is 
that the effective coupling between the sky emission and 
the inner error beam (^ 5') of the Mopra telescope is 
not as uniform as is implicitly assumed by the use of r]c 


= 0.55-0.57 in Table The conversion suggests instead 
that an efficiency of 0.50 may be better value, i.e., about 
10% below 77c, but 15% above 775 . This would be true 
if the sky emission only partly fills (perhaps to ^60% 
of) the inner error beam, giving an effective 77 between 
these nominal values. However, the appearance of the 
maps (especially the from which this result is ob¬ 

tained) is not completely consistent with this idea: over 
large areas the cloud structure is very extended, al¬ 

though it is compact in other areas, so we might expect 
a more variable conversion between the ThrUMMS and 
CfA brightness scales, rather than the fairly systematic 
trend evident in Figure 

An additional factor may be that our BSM technique 
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Figure 10. Mosaic of DR3 integrated intensities (zeroth moment) for the three CO isotopologues across ^^9° of longitude within the 
ThrUMMS survey area. Here, red = ^^CO, green = ^^CO, and blue = C^^O. Note the wide variation in line ratios across the Galactic 
Plane, indicated by the varying colors in this image. E.g., the prominent red cloud at I ~351°5 is relatively bright in ^^CO, but relatively 
faint in and C^^O. In contrast, there are numerous small “blue” features throughout the image, where is relatively bright, but 

and 12CO less so. 


is systematically missing a few percent of the line flux 
from all mapped areas, due to the undersampling of the 
emission structure (see 12.1.2). That is, between raster 
lines there may be small, bright, compact (i.e., beam¬ 
sized) structures whose total flux would not be fully 
measured by the beam-sampling, yielding an apparently 
lower brightness scale for the maps after interpolating 
between the raster lines. Another way to understand 
this issue is by analogy with interferometry: there are 
scales to which our maps are less sensitive than if they 
were properly Nyquist-sampled. Tests of this effect in 
maps of the CHaMP areas (which are Nyquist- 
sampled; Barnes et al, in prep.) suggest that it may be 
relevant, but we defer that discussion to that work. For 
the purposes of this paper, we use these explanations as 
working assumptions to account for the 0.90 conversion 
factor. Thus, while we quote all results herein on the T* 
scale as derived with the tabulated r]c factors, the reader 
should be aware that these brightness levels may need to 
be increased by a further ^10%. 


3. DATAPRODUCTS 
3.1. Online Aeeess of Speetral Line Cubes 

From its inception, ThrUMMS was designed to be 
an open-access project, with frequent Data Releases as 
progress is made in the observing and data reduction. As 
of 2014 May, we have placed a large amount of data on 
the website http://www.astro.ufl.edu/thrumms, freely 
available to the community. For the and 

cubes, the data are at DR3 level and organised into 6°- 
wide “Sectors” (i.e., each including up to 48 contiguous 
half-degree fields; see Figs. [^11). For the and CN 

data, we currently have avanaBTe the DR2 cubes, organ¬ 
ised into 2°-wide blocks; DR4 is planned to include the 
augmented artifact-correction for the latter two species, 
as well as wider area coverage from the 2014 observing 
season, and is anticipated to be available in mid-to-late 
2015. At the present time, these data cubes cover ^80% 
of the 60° X 2° target area for all 4 species. With the 
anticipated completion of observing in 2015, later Data 
Releases will continue to be provided as processing con¬ 
tinues and further data products can be assembled. 


3.2. Moment Maps and Initial Proeessing 

We use the Miriad package (jSault et al.|1995|) to pro¬ 
cess the data cubes into moment maps and LV diagrams, 
and perform other basic analyses, with visualis ation pro- 


4ween 41 


Figure 
hrUMMS 


vided by the Karma package (Good 
gives some sample comparisons between 
maps and other surveys, while Figure shows sample 
VFM data for our 4 molecular transitions. Note the high 
S/N of the and maps, and even though the 

and CN maps have noticeably smaller S/N, the 
peaks in these maps (around the positions of the bright¬ 
est or emission) are still detected with good 

reliability (S/N > 5-10). In Figure we show further 
examples of higher moment and LV maps in the 
line, demonstrating the very high image fidelity achieved 
across extended emission regions. 


3.3. DR3 and Smooth and Mask 

The example moments maps of Figurej^were processed 
with a smooth-and-mask (SAM) algorithm in order to 
improve the fidelity of these products. As part of this 
algorithm, the DR3 cubes are binned by 4 channels in 
velocity (see Table |^, which improves the S/N in the 
public data cubes and reduces the file sizes. (The full- 
resolution cubes are, however, still available on req uest.) 
Our SAM al g orithm is simi l ar to t hos e described by [Rots 


et al. (1990 

), Dame et al. ( 

12001 

), or |Dame (2011) 


form a data cube that has been smoothed/convolved in 
all three dimensions by a factor of 2 (i.e., beyond the 
4-channel binning of the public data), a 5crrms thresh¬ 
old was calculated, and voxels (where a voxel is defined 
here as a single (/,5,U) element) below this level in the 
smoothed cube were used to define a blanking mask for 
voxels in the unsmoothed data cube. Moments are then 
calculated only with the unmasked voxels. 

Construction of moment maps, however, depends 
strongly on the particular investigation one wishes to 
pursue with the data cubes. There is a large range of 
input parameter combinations possible for such maps, in¬ 
cluding angular extent, velocity range, projections, SAM 
parameters, uncertainty estimation, etc. For this reason, 
moment maps using the above techniques and based on 
the public data cubes are not currently part of the DRs, 
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Figure 11. Mosaic of ThrUMMS pre-DR 4 l- V images (SAM-based zeroth moment integrated across all b) for (red), (green), 

and (blue) across the 4Q. As in Fig. |lO[ we see wide variations in line ratios across the Galactic Plane from the varying colours in 

this image. The data maximum, saturated colour, and 3a (almost black) levels in this image are respectively at 577, 248, and 12Karcmin 
(^^CO), 187, 63, and 5Karcmin (^^CO), and 37, 5, and 3Karcmin (C^^O). 


but are available upon request (such requests are pro¬ 
cessed on a best-efforts basis). 

We do plan, however, to include sets of standardised 
moment maps in future DRs, once we have re-observed 
most of the weather-damaged areas and made the data 
cubes’ noise statistics more uniform. Meanwhile, users 
can use many standard astronomical or image processing 
packages to perform their own analyses, including com- 



The full survey data cubes provide an extremely rich 
dataset for current and future investigations. This rich¬ 
ness means it is difficult to render, whether on the printed 
page or a web browser, the detailed physical insights 
contained in the data. Nevertheless, the reader may 
gain some intuitive feeling for this richness by means of 
colour-compositing the different species’ integrated in¬ 
tensity maps. We show examples of this in Figure [Tol 
combining the DRS ^^CO, ^^CO, and integrated 

intensities, and in Figure pT| for pre-DR4 versions of the 
^^CO, ^^CO, and data, each showing clear posi¬ 

tional variations in the line ratios as variations in the 
rendered colours. 

Similarly, compressed image versions of the spectral 
line integrated intensities across the full survey area, in 
all four species, are given in the Appendix. The reader 
should note, however, that the Appendix images, and 
the data used in the analysis presented below, have been 
manually edited from the DRS versions available on the 
ThrUMMS web site to exclude most of the artifacts 
generated during poorer observing conditions (“weather 
stripes”). These excisions are manifest in the Appendix 
figures as half-degree-wide gaps in certain fields; the gaps 
themselves are not the weather stripes. These blanked 
areas cover only ^2% of the mapped area, however. 

4. FIRST SCIENCE RESULTS 

ThrUMMS is designed to provide legacy data for a 
wide array of science investigations, with specifically- 
focused applications in astrochemistry and the physical 


state of GMCs; the neutral/molecular ISM connection; 
assessment of the turbulence and dynamical parameter 
space on both a targeted (i.e., specific regions) and holis¬ 
tic level; and Galactic structure and kinematics. But 
the quality and scope of the dataset is also suitable for 
studies where our attention is not currently focused. Be¬ 
cause ThrUMMS is an open project, we welcome new 
initiatives from members of the community not yet in¬ 
volved with ThrUMMS, and indeed encourage those who 
are so inclined to collaborate with us in developing new 
software or data products to expand ThrUMMS’ impact 
beyond our original vision. 

Our first detailed analysis of a specific region, including 
derivation of kinematic distances, the holistic properties 
of the GMC population, their relationship to the HI and 
cold dust (Hi-GAL+MIPSGAL) emission, and implica¬ 
tions for the complexes’ inferre d star formation rat e, is 
given in Paper II of this series (Nguyen et al.||2014) for 
the G333/RCW106 complex. In contrast, we describe 
here some first results of the global analysis of the CO- 
isotopologue lines across the 4Q. 

4.1. Optical Depth and Excitation Temperature 
Variations 

With our comprehensive coverage of a large portion 
of the Milky Way in multiple species, we have an un¬ 
precedented opportunity to use the molecular lines to 
investigate global GMC conditions. We examine now 
the variation in the CO line ratios (and consequently, 
the variation in the physical conditions in GMCs) with 
position in the Galaxy, as shown visually in Eigures pR| 
and 11 This rendering can be intuitively understood in 
terms of radiative transfer physics in the following way. 

The clouds rendered in redder colours show where 
the ^^CO is relatively much brighter on average than 
either ^^CO or C^^O, i.e., where the ^^CO/^^CO and 
^^CO/C^^O brightness ratios are near their maximal val¬ 
ues. This must therefore indicate areas of relatively low 
optical depth, since the intrinsic abundance ratios are 
also quite large (any line saturation = high optical depth 
would produce brightness ratios significantly less than 
the respective abundance ratios). At the same time, the 
lines in such areas are quite bright, typically with 
T* > 30K. So the “red” clouds are typically those of low 
r and high Tex, or “warm and translucent” clouds. 
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Figure 12. (left) Plot of voxel density (coded by color and contoured on an arbi trary scale) in the Ii3(jq-Ii2(jq plane, where the intensity 
is integrated over lkms“^ velocity bins, from all-4Q mosaics similar to Fig. |lO| The lines are labelled by ri 3 for an assumed Ris = 60. 
(right) Similar plot for and ^^CO. The lines are labelled by rig for an assumed Rig = 10. See text for further discussion. 


In contrast, green or even cyan/blue areas are where 
either ^^CO, C^^O, or both are much brighter compared 
to than on average. This is only likely to occur 

where the optical depth is high enough to produce some 
level of saturation in either line, so that the ^^CO/^^CO 
and ^^CO/C^^O brightness ratios are closer to unity than 
average. At the same time, the emission in such 

areas is not particularly bright; coupled with the higher 
optical depth, this means that the excitation tempera¬ 
ture in such clouds is likely to be somewhat low. So the 
“green” and “blue” clouds will typically represent areas 
of high r and low Tex, or “cold and opaque” clouds. 

Therefore, we understand that the line ratios provide 
a simple yet powerful way to derive basic physical con¬ 
ditions in Galactic GMCs from an elementary radiative 
transfer analysis. The rendered colours make these re¬ 
lationships intuitively clear: they are set so “white” or 
“neutral” colours represent the average line ratios and 
not where the line ratios are close to 1. 

Other workers have used a similar approach to derive 
cloud properties from GO line rat i os; ex ample s include 


Myers et al. (1983), Wong et al. (2008), and Hindson 
et al. " pOT? ). Eacn of these, however, uses some sim- 


plifying assumptions about the excitation temperature, 
abundance ratios, or other radiative transfer parameters, 
and typically for relatively small samples of clouds cov¬ 
ering at most 10s of pc. Ours is the first such study of 
which we are aware that attempts this analysis uniformly 
across many kpc of the Galaxy; at the same time, we re¬ 
tain the pc-scale resolution necessary to characterise the 
clumps that give rise to individual star clusters. 

In the analysis that follows, we use the full 3D iso-GO 
data cubes, but binned in velocity to lkms“^ channels, 
to compute all physical quantities. We do this in order 
to fully track the variations of Tex- r, and other parame¬ 
ters with position and velocity, but also to boost the S/N 
in the C^^O data especially. In what follows, therefore, 
the notation JdV refers only to these lkms“^ integrals, 
which are still voxels in the binned cubes, unless other¬ 


wise stated. 

For abundance ratios Ris = [^^CO]/[^^CO] and Rig 
= [^^CO]/[C^^O], the optical depths in the different iso- 
topologues will follow ri 2 = RisTis and ris = RigTis- 
Then the radiative transfer equation for each spectral 
line (in the Rayleigh-Jeans limit, where Jj^ cx: T), 


Fmb = (^ex - ^bg)(l - e 


( 2 ) 


can be combined to convert the line ratios to an optical 
depth, with the further assumption of a common exci¬ 
tation temperature Tex- In general, this last condition 
may not hold everywhere: the lower optical depth 

means that this species will more often be subthermally 
excited compared to ^^GO, which has a low er “effective' 


density due to stron g radiative trapping (Evans 1999, 
Heyer fc Dame|2Q15 ). Nevertheless, this assumption wil 


only act to unde restimate the resultant column 

density (see ^4.2). 

With these assumptions, the brightness ratio in the 
ThrUMMS maps should obey 


^ 1 - 

Ti2 1 - ’ 


(3) 


and a similar equation applies for and ^^CO. For 

typical values of R and r (especially for Ris), the de¬ 
nominator in eq. (3) will often be close to unity. Also, 
for any reasonable value of Ris, the will be virtu¬ 

ally opaque almost everywhere, and T 12 = (Tex — Fbg). 

The visual representations of line ratio variations, as 
seen in Figures pT] a nd 0 can be quantified and under¬ 
stood in FignreTT^ Here we see that, in the left panel, 
there is a very wide range of ratios log(/i3co/^i2co) = 
log(/ TiscpdV/ f Ti 2 codV), from ^ +0.2 in the upper- 
left of this panel (seen as green areas in Fig s.[T^ and 
pr] ) to ^ -1.1 in the lower-right of Figure fl^fred- 
areas in Figs.j^ and 11). Similarly in tHe right 
panel, we see a range of ratios for log(/ci8o/T3co) = 
log(/ Tc isodF/ f Ti3QodF), from +0.5 in the upper- 
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Figure 13. Two plots of the same CO line ratio data, {left) Voxel density (coded by colour and contoured on an arbitrary scale) in the 
CO line ratio plane. Overlaid are simple radiative transfer models based on eqs. (2-3); see text for further discussion, (right) The same 
CO line ratio plane as in the left panel, but the colours in this plot represent the approximate appearance of pixels in composite images 
like Figs.|10| |11| and in the Appendix. The overlays are the same as in the left panel. 


left of this panel (seen as blue or cyan areas in Fig. 10) 
to ^ -1.2 in the lower-right of Figure [T^ (green or yellow 
areas in Fig.fW 


Now we can apply eq. (3) to the data in Figure 12 
where in each panel we show lines of constant r for an 
assumed R. As expected, low optical depth for the line 
pairs occurs to the right of each panel, while high optical 
depth lies to the left. Interestingly, in Figure [T^ we see 
several groups of high-brightness emission peaks(i.e., the 
brightest clumps in Figs. 10, 11) as separate extensions of 
points, parallel to lines oiconstant r, in the distribution 
to the upper right of each panel. These extensions appar¬ 
ently correspond to groups of bright cl ump s (complexes) 
with different regional r and Tex (see ^4.3[ ). 

The plots in Figure can be cornbined into a sin¬ 
gle “CO-colour” plot (Fig. 13), where a voxel’s position 
in this plot corresponds to its colour rendering in images 
like Figure or This can be compared with the ana¬ 
lytical relations derived above (eqs. 2-3) which are shown 
as overlaid lines in Figure From this comparison it 
is immediately obvious that we cannot model the most 
extreme “green” and “blue” voxels with the assumptions 
used above (i.e., a common Tex and fixed ratios R>1). 
For this small number of voxels, the model would require 
either (1) ratios R<1 (e.g., due to strong fractionation 
or other selective effects), or (2) a lower Tex in the more 
abundant and optically thick species (e.g., due to a colder 
foreground layer and warmer cloud centre). Fractiona¬ 
tion effects for both and especially (|Glassgold 

et al.||1985) do occur, as does self-absorpton from colder 
cloud envelopes, so undoubtedly these both contribute 
to deviations from the simple radiative transfer picture 
presented here. Note also that single clouds may be rep¬ 
resented by a range of loci in Figure So, while we 
can understand most of the line ratio variations in Fig¬ 
ures!^ and i.e., the lower and leftmost boundaries of 
the distribution in Figure eventually a more sophisti¬ 
cated radiative transfer approach will be needed to fully 
understand all the line ratios, because the interplay be¬ 
tween optical depth, excitation, and abundance is more 


complex than the above treatment allows. 

Nevertheless, we can use our treatment to estimate the 
molecul ar c olumn density distribution in the Galactic 
Plane ( ^4.2[ ). We use the line ratios to evaluate r via 
eq. (3) at each 1 kms“^-binned voxel, and then invert 
eq. (2) to also obtain Tex at each voxel, producing all- 
4Q cubes of both parameters. Finally, we form moment 
maps from each of these cubes in the normal way, e.g., 
integrating over large velocity ranges to obtain the total 
r, peak or average Tex, or total Nqo (see also ^4.2). We 
show examples of these results in Figure although 
these are moment maps, the structures they exhibit are 
typical of the cubes over which they were integrated, al¬ 
though there are fewer visible complexities in any given 
velocity channel. 

The r panel in Figure (top) shows a very high- 
contrast image: all the highest optical depth structures 
are quite compact, typically having size scales of ^1'- 
5'. It may not be coincidental that this is also the size 
scale (^a few pc at 3kpc distances) of dense clumps, 
IRDGs, and clusters. Lower optical depth structures are 
not quite so spatially compact, and are somewhat widely- 
distributed. 

In contrast, the Tex panel of Figure (middle) shows 
a dramatic difference: nearly all structures are rela¬ 
tively large-scale, and somewhat amorphous and/or low- 
contrast. While the warmest structures are quite ex¬ 
tended, at least 5', moderate-temperature clouds extend 
over even larger degree-wide (^50 pc at 3kpc) scales. 
Again, it is tempting to draw a connection between these 
scales and the similar sizes of GMCs and classic HII re¬ 
gions, 10-100 pc. 

We note an interesting consequence of this radiative 
transfer approach. In Figure rig maximises near a 
value of 0.6; for any higher optical depth, all GO iso- 
topologues would begin to saturate at T^b = ^ex Fbg. 
Such clouds would have r indistinguishable from infinity, 
whatever the Tex (all of the Rig curves converge to where 
both line ratios are 1). Even so, rig = 0.6 corresponds ap¬ 
proximately to ri 2 ^ 250 for any reasonable values of the 
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335 ® 334 ® 333 ® 332 ® 331 ® 

Goloctic Longitude 


Figure 14. Sample integrated maps of ri 2 (top, purple), Tex (middle, orange), and Nqo (bottom, blue) near I = 333°, obtained from the 
respective cubes of these quantities as described in the text. Note in particular that t he c olumn density map is noticeably more filamentary 
than either the very clumpy r map, the rather amorphous Tex map, or the Ico (Fig. |18| > map. 


abundance ratios R. This means that there are a large 
number of clouds in the ThrUMMS maps with very large 
column densities (e.g., near the top-right of the plots in 
Fig.f^ with blue, cyan, or green ratios), some of which 
could^asquerade as lower column densiy clouds if one 
were to convert from /12 with a standard X-factor. 

4.2. Column Density Maps 


A subject of much study in the literature is the con¬ 
version of observed molecular emission into estimates of 
molecular column density or mass, which are important 
for a number of questions in astrophysics, e.g., in star 
formation when measuring cloud stability against gravi¬ 
tational collapse, or across disk galaxies for studies of the 
Kennicutt-Schmidt relation, or for estimating the contri¬ 
bution of the molecular ISM to a galaxy’s dynamics, etc. 


Tex [K] 
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Dame et al. (2001) summarised this topic and, when con- 
verting CO emission into mass, obtained a best average 
estimate for the X-factor of 

Xco = Nh2 / Ico 

= 1.8 X 10^^ H 2 molecules m“^/(K km s“^) (4) 

in the solar neighbourhoo d, based on high-la titude data 
in the CfA survey. But asjDame et al.|(|2001|) point out, 
this value is expected to vary locally in the Galactic mid¬ 
plane (i.e., at low latitudes) due to various factors. 

Having derived r and Tex at each voxel, we can use the 
ThrUMMS data to make a unique contribution to this 
exercise with a CO column density map, via 


Nqo = 


3h Q(Tex)e^“/'=^» 


Jt,(l — 


/ 


rdU, (5) 


where /i is the CO dipole moment, Q is the rotational 
partition function, Eu and Ju are the energy and quan¬ 
tum number of the upper level of the transition at fre¬ 
quency and the integral is over the velocity range of 
the emission line. Here, instead of using a constant Xqo 
factor to convert the observed Iqo into a molecu lar gas 
mass, as is so metimes done in the literature despite [Dame 
et al. (2001 )’s caveat ag ainst doing so, we can use the 
(r,'iex) estimates above (^4.1) to compute a more realis¬ 
tic CO column density. Tnis can then be compared with 
other estimates of molecular gas mass (e.g., SED fits to 
the Hi-GAL data) to derive CO abundance maps. 

The computation of eq. (5) across a spatially-resolved 
map with ^8x10^ independent (/,5,U) resolution ele¬ 
ments, complete with a calculation of the partition func¬ 
tion at each voxel and appropriately propagating uncer¬ 
tainties in each quantity, is a non-trivial exercise. More¬ 
over, in many locations the ^^CO and especially C^^O 
can be quite weak, so there are S/N limitations in apply¬ 
ing this method to our full -reso lution data cubes. Nev¬ 
ertheless, as described in §4.1| we have performed this 
computation on cubes that have been binned in veloc¬ 
ity to lkms“^ channels. After so computing the r and 
Tex in each such voxel, we use eq. (5) to derive cubes of 
column density per channel. 

We also describe here a simpler calculation, which 
gives a very good approximation to eq. (5) in cases 
where insufficient computational power is available. The 
temperature-dependent terms in eq. (5) can be collected 
and written as a closed function of Tex, which we call the 
modified temperature function Tmod- 


^mod — 


Q(^ex)' 




Ju{^ — e ku/kTe^^ 


(6) 


Over a typical range of Tex = 5-100 K, the following for¬ 
mula for Tmod gives deviations of <1% from the exact 
expression (eq. 6), without having to calculate the parti¬ 
tion function Q explicitly: 


log"(2:„,od) = 1 + log" [0.053 . 


( 7 ) 


Having computed the Tmod cube from the Tex cube 
using either eq. (6) or (7), we then obtain the column 
density cube from this and the r cube. 


Nco = 5.57 X 10^® m ^ Tmod / rdV , 




(8) 


where the integral over V means per lkms“^ channel. 
We emphasise, however, that we have actually done the 
full calculation from eq. (5/6) in this study. 

A moment map of the column density cube from eq. (5) 


is shown in Figure 14 (bottom panel). This no w embod- 


ies a muc h larger ^Tonversion factor” than in Dame et] 


al. (2001). In the cool and opaque case, we might have 
typical values r>100 and Tex^lOK. For a cloud with an 
effective line width ^5 kms“^, Nqo ^ 4.1x10^^ m“^. In 
fact, the bottom panel of Figure 14 shows that column 
densities computed in this way can reach a peak value 
some 5x higher than this, or Nco ^ 2xl0^^m“^. As¬ 
suming that the abundance is 10“^ (although see 

^4.3), this is equivalent to ^ 0.4-2x 10^^ m“^, or a 
mass column F ^ 760-3700 Mq pc“^. Even the lower end 
of this range is close to the mass columns in p arsec-scale 
“d ense c lumps” seen in previous surveys (e.g.. 


Barnes et 


ah] 2011), but is much higher than typical GlVlC mass 


columns. We can convert the above cool, opaque exam¬ 
ple to an equivalent Xco^ 8-40x10^^ m“^/(Kkms“^), 
and so realise that the Milky Way’s molecular gas mass 
may need to be revised upwards by a significant amo unt, 
at least in some locations; see the discussion in §4.3| 

Equally remarkable as the actual numbers are the 
structural revelations in Figure As has been seen 
in many previous studies, the integrated intensity iso-GO 
maps (Appendix, Fig.[TM^ have a rather complex com¬ 
bination of structures. I'he respective data cubes reveal 
that this complexity breaks down into individual regions 
or clouds at various Wsr, which may be clumpy and 
compact, or filamentary, or amorphous on larger angular 
scales. However, the cubes by themselves do not reveal 
any obvious systematic trend among these structures. 

In contrast, the r, Tex, and Nco cubes show system¬ 
atically d iffer ent structural characteristics. As already 
noted in §4.1[ the r cubes are clumpy and high-contrast, 
while the Tex cubes are amorphous and low-contrast. In 
Figure 14, the A^co panel (bottom) reveals the combi¬ 
nation oTthe r and T^xi a high-contrast and more fila¬ 
mentary medium than either of the above constituents. 
One is naturally drawn to compare such images with the 
almost ubiquitous filamentary structures evident in the 
Hi-GAL maps. 


4.3. Global Implications 

Based on the distribution of points in Figure and 
the conversion embodied in eq. (5), we can in principle 
compute mass, density, column density, or any other dis¬ 
tribution functions for the entire molecular cloud popula¬ 
tion of the 4Q, such as in Figure [T^ (discussed next). For 
many derived parameters (such as mass), we also need 
distances to individual clouds to obtain them from the 
column density computed here. However, assignment of 
distances to all the structures we see in the cubes is far 
beyond the scope of this paper. 

Despite this, our results have a number of significant 
implications for studies of molecular cloud populations 
in disk galaxies. For example, a variety of observa¬ 
tional measurements and chemical models indicate that 
the fractional CO abundance relative to molecular hy¬ 
drogen typicall}!^ lies in the range 10“^‘^ to a 

This is an average across different environments throughout 
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10 lO'lO^O^O^O^ 

li2co [K.kms“^] Incidence 


Figure 15. Voxel incidence of the ratio of computed column density from eq. (5) to the measured integrated intensity per lkms“^ 

binned channel, plotted as a function of the integrated intensity, across the 4Q. The green line represents mean values of this ratio in 
narrow bins of Iqo] the dashed blue line shows the expected Xco/Ri2 ratio, or equivalently, the conversion of Ico to using eq. (4) 
with standard values of X and R] while the sloping blue line is a weighted power-law fit to the modal trend, as labelled at the top (see 
text for further discussion). The histogram to the right is effectively a distribution function for XjR: a column-weighted mean value is 
5.1x10^° CO moleculesm“^/(K kms“^). 


mean value of 10“^ is often used in the literature. As 
discussed in the last section, this means that we should 
be able to simply convert our computed Nqo to A^H2 by 
such a factor. But we also have the standard conversion 
of Ico to A^H 2 using eq. (4). For these two conversions to 
be consistent, the various quantities must be related. To 
make these relationships clearer, we rewrite eq. (4) as 

A^H2 = -^CO^CO = ^12^^00 5 (9) 

where we define the ratio R 12 from the actual CO 
abundance in the gas, R 12 = A^H 2 /^cO: or typically 
2Q+3.5to+4.5 above. Note that this definition makes 
the gas abundance ratio R 12 different to the empirically 
determined Xqo conversion factor. Then 

Xco/Ico =-^co/^i2 • (10) 

With the standard values on the RHS of eq. (10), the 

the Milky Way: lower values seem appropriate in dark 

clouds — our cool, opaque case — and higher values in lu¬ 

minous HII regions — our warm, translucent case (e.g., see|Gerner| 
|et al.|2014| >. - 


ratio on the LHS should be a constant, 

Xco/Ico = 1-8 X 10^^ CO molecules m“^/(K kms“^). 

( 11 ) 

However, when we plot eq. (11) vs. Ico iu Figure [l^ we 
see a marked variation in the actual ratio X/R as shown 
in eq. (10), instead of the flat line that would be expected 
from eq. (ll) with fixed values of X and R (shown as a 
dashed blue line in Fig.p!^. 

Even at the lowest values of Ico (~6Kkms“^ per 
lkms“^-wide channel) in Figure [T^ the modal fit (blue 
line) of X/R already begins above the standard value 
(dashed line), while the mean values (green line) lie 
around 3-4x10^^ CO molecules m“^/(K kms“^), or 
twice the standard value. If one assumes Xco is fixed 
to the standard value, this implies that R 12 = 5x10^, 
or half its standard value. On the other hand, if one 
assumes R 12 is fixed to its standard value, this implies 
that Xco = 3-4x10^^ H 2 molecules m“^/(Kkms“^), or 
twice its standard value. As Ico rises, the modal fit 
hi X/R also rises further to almost 10^^ CO molecules 
m“^/(Kkms“^) at Ico = 40Kkms“^ per 1 kms“^-wide 
channel. This means that at these brightness levels, ei¬ 
ther Xco is 4x higher than standard, or Rco is 4x lower 
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than standard, or some combination of these two options 
applies. At Iqo > 20 kms“^, the Xco /^12 distribution 
actually bifurcate^^into two regimes, one with the ratio 
up to 2x higher than the standard value, and another 
up to 10 X higher. A weighted power-law fit to the upper 
trend (sloping blue line in the plot) gives 

X/R (X (12) 

across the measured range of /co? or equivalently, 

Nco = 1.6 X 10 "° m-2 /1-38. (13) 


This is a remarkable result for several reasons, but 
mainly because in retrospect, it should have been more 
widely recognised before now. For example, the fact that 
eq. (13) is superlinear falls naturally out of our very sim¬ 
ple radiative transfer ass umptions, and theoretically was 


predicted on this basi s (Krumholz & Thompson 


2007 


Narayanan et al.|[2008 ). Also, because eq. (13) is super 


linear, it is no longer as physically meaningful to talk 
about a single X-factor: the direct conversion from CO 
integrated intensity to hydrogen column density can be 
seen as an illusion that is correct only in the mean, and 
masks the very real effects of high optical depth in the 
line. Thus, the scale factor in eq. (13) has units 
of molecules m“^/(K kms“^)^‘^^, and so dimensionally 
is not a “conversion factor.” 

Nevertheless, we can compute average quantities 
from our data, which may still be useful in non- 
spatially-resolved studies, such as extragalactic or high- 
redshift star-formation settings. The equivalent me¬ 
dian value of the X-factor for all our data is Xqo = 
3.8x10^^ m“^/(K kms“^) or 2.1 x the standard value, 
while a column-weighted median conversion is even larger 
than this, Xqo = 5.1x10^^ m“^/(Kkms“^) or 2.8x the 
standard value. 

However, eqs. (12-13) represent just the modal trend 
in Figure [l^ It should be clear that the range of the 
variations in Xco /^12 is far too large to be attributed 
to changes in just one of these parameters. Even the la 
variation from the modal trend in Figure [Ts] is 0.3 in the 
^ ~ ~ ' T5]show 

ler than 


log, or a factor of 2. But the extrema in Figure 
X/R ratios from about 30 x larger to 100 x sma 
standard. This represents >3 orders of magnitude em¬ 
pirical uncertainty in the conversion of Iqo to . This 
variation exceeds any uncertainty in the assumptions we 
have used in this analysis, such as the assumption of LTE 
at a single Tex, or values of the ratios T, or even of the 
constancy of Xqo in the face of the metal licity gradient 
in the Galactic disk (|Heyer & Dame||2015). 

Such uncertainties aside, in the mean it is likely that, 
as Iqo rises from 4 to 40Kkms“^ in each lkms“^-wide 
channel, Xqo is very likely rising and that R 12 is prob¬ 
ably falling (meaning that the gas phase abundance of 
relative to H 2 is probably rising). Therefore, the 
conversion of Iqo to an equivalent mass surface density is 


T his bifurcation, which also appears in the left panel of 
Fig. |12| appears to be entirely due to the two brightest GMC com¬ 
plexes in in our ma ps, b ut these still illustrate how the dis¬ 

tribution of points in Fig.[^ depends on the iso-CO ratios. The 
upper branch comes from the G333 complex, which is bright in all 
three GO lines, and so has moderate optical depth and high column 
density. The lower branch comes from the G351.3+0.6 complex, 
which is bright in ^^GO only, and so has lower optical depth and 
lower column density. 


not as simple as has often been assumed. While we can¬ 
not yet derive the exact corrections to cloud masses, and 
the Milky Way’s total molecular mass, without clouds’ 
individual distances, this result has obvious implications 
for studies of the molecular mass distribution in the 
Milky Way, star formation laws in disk galaxies s uch as 


the Schmidt-Kennicutt relations, and other topics (Leroy 
et al.|2QQ8 Lada|2Ql^ ). At the very least, it implies that 


such relations must be recalibrated for the highest optical 
depth clouds and regions, which in our data are widely 
distributed across the Galactic Plane, and seem to have 
not yet been adequately taken into account. 


5. GONCLUSIONS 

We have presented information about, and first data 
releases and science results from, a comprehensive new 
mm-wave multi-spectral-line survey of the southern 
Milky Way, ThrUMMS. Essential features of this sur¬ 
vey are its public access, speed, coverage (120 deg^), 
high angular (arcminute) and physical (parsec) resolu¬ 
tion, and simultaneous mapping of four important molec¬ 
ular tracers, ^^CO, ^^CO, G^^O, and GN. This combina¬ 
tion of attributes makes ThrUMMS unique. We have 
also presented first science results on the global prop¬ 
erties of the Galactic population of molecular clouds in 
the Eourth Quadrant. These include the discovery of 
widespread and large variations in the line ratios of the 
three main GO isotopologues with position and velocity, 
and the implication that some of these variations arise 
in a large population of clouds with very high column 
density. We derive a new calibration of the CO-to-H 2 
conversion, embodied by our eq. (13) and replacing the 
concept of a fixed Xco- With an eventual catalogue of 
cloud distances, this will drive related corrections to star- 
formation and other laws. We anticipate a wide range of 
further science applications to spring from this unique 
data set for many years into the future. 
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APPENDIX 

THRUMMS INTEGRATED INTENSITY IMAGES 

We present here all integrated intensity images across the 4Q survey area, broken up into one 10° x 2° sector per 
printed page. Each sector is shown both as a combined pseudo-colour image with red = ^^CO, green = ^^CO, and 
blue = C^^O, and beneath this, four separate images for each of the four species mapped, as labelled. A single 
contour at 28.94 K kms“^ ~ 20cr is shown in each of the four species’ panels, labelled in the panel. Note that 

CN (yellow) is not included in the pseudo-colour images. 
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Figure 16. ThrUMMS DR3 integrated intensity images for 360° > I >350°. 
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Figure 17. ThrUMMS DR3 integrated intensity images for 350° > I >340°. 
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Galactic Longitude 


Figure 18. ThrUMMS DR3 integrated intensity images for 340° > I >330°. 
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Figure 19. ThrUMMS DR3 integrated intensity images for 330° > I >320°. 
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Galactic Longitude 


Figure 20. ThrUMMS DR3 integrated intensity images for 320° > I >310°. 
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Figure 21. ThrUMMS DR3 integrated intensity images for 310° > I >300°. 





















